Vibrational edge modes in intrinsically heterogeneous doped transition metal oxides 
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By applying an unrestricted Hartree-Fock and a Random Phase approximations to a multiband 
Peierls-Hubbard Hamiltonian, we study the phonon mode structure in models of transition metal 
oxides in the presence of intrinsic nanoscale inhomogeneities induced by hole doping. We identify 
low frequency local vibrational modes pinned to the sharp interfaces between regions of distinct 
electronic structure (doped and undoped) and separated in frequency from the band of extended 
phonons. A characteristic of these "edge" modes is that their energy is essentially insensitive to 
the doping level. We discuss the experimental manifestations of these modes in inelastic neutron 
scattering, and also in spin and charge excitation spectra. 



Recent advances in experimental techniques provide 
strong indications that a broad class of electronic ma- 
terials, commonly referred to as strongly correlated, eyi- 
hibit intrinsically heterogeneous electronic phases^, 0, 
Hi Q IE IE HI' Strong interactions in these materi- 
als drive phase separation between regions of distinct 
electronic structure 9]. However, a competing interac- 
tion or electronic kinetic energy can lead to frustra- 
tion of the global phase separation and to formation of 
nanoscale inhomogeneities, which can take the form of, 
e.g., stripes/checkerboards 10] or clumps|lj|. Unlike the 
conventional spin and charge density wave phases, the or- 
dering of intrinsic heterogeneities is sensitive to thermal 
and quantum fluctuations T^l, as well as static disorder. 
This makes an unambiguous identification of such phases, 
and their distinction from the more familiar density wave 
instabilities problematic, requiring development of new 
diagnostics more sensitive to the unique structure of such 
phases. Here, we propose one such diagnostic tool which 
relies on detection of localized edge modes pinned to the 
interface between the regions of different electronic struc- 
ture. The edge (interface) modes appear in every de- 
gree of freedom coupled to the inhomogeneity (charge, 
spin, lattice) jl^ . with corresponding correlated experi- 
mental signatures. Here we study the lattice vibrational 
edge modes which lend themselves to study with inelastic 
neutron scattering. We also compute corresponding ex- 
citation spectra in spin and charge channels and describe 
the cross-correlations. 

We focus on one particular class of transition metal 
oxides that has attracted significant attention recently, 
namely cuprates. These are the materials exhibit- 
ing high-Tc superconductivity upon doping. The pres- 
ence of stripes/checkerboards and their possible influ- 
ence on superconductivity are subjects of intensive de- 
bate. Since it appears that stripes can significantly af- 



fect suDerconductivitv [l3. Iisl Il6|. it is important to es- 
tablish whether they are indeed present in the materials. 
Our purpose here is to elucidate typical manifestations of 
stripe order in experimental probes sensitive to the exci- 
tations in spin, lattice and charge sectors. In particular, 
we compare our theoretical results with the observations 
of "anomalous" phonon modes in the cuprates [itI fl^ . 
Analogous results also apply to edge mode signatures in 
other doped transition metal oxides (and related com- 
plex electronic materials), as found in recent experimen- 
tal studies [iolE^ Esl l24l| : corresponding modes (phonon 
shape modes) have already been predicted and identified 
in conjugated polymers and quasi-lD (one-dimensional) 
charge transfer solids [IE 113 . 

To model the Cu02 planes of doped cuprates, we 
use a 2D three band extended Peierls-Hubbard Hamilto- 
nian, which includes both electron-electron and electron- 
phonon interactions [IE HE]: 



Ho= tpdiuij)ictcja + H.c.) + J2'^^i^ij)4a<^^<y 

<ij>(7 i,fT 
<ij> i,j,{T,a' 

Here, creates a hole with spin cr on site i; each 
site has one orbital {(1x2 _y2 on Cu, or O p^, or pj, on 
O). The Cu (O) site electronic energy is (cp). Here, 
Uij represents the on-site Cu (O) Coulomb repulsion, 
Ud {Up), or the intersite one, Upd', and the summation 
J2 is taken except for the case («,cr) = {j,<^')- The 
electron-lattice interaction causes modification of the Cu- 
O hopping strength through the oxygen displacement Uij : 
tpd{uij) = tpd ± auij, where ) applies if the Cu-0 
bond shrinks (stretches) for a positive Uij; it also af- 
fects the Cu on-site energies f-d{uij) — ed + (i^uOi 
where the sum runs over the four neighboring O ions. 
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The other oxygen modes couple to electron charge more 
weakly, and are neglected for siniplicity. We use the fol- 
lowing set of model parameters |la l27l |: Cp — = 4.4 
eV, Ud = 11 eV, Up = 3.3 eV, Upd = 1-1 eV, and 
K = 38.7 cV/A2, a = 5.2 eV/A, (3 = 1.2 eV/A, with 
tpd = 1-1 cV. To approximately solve the model, we 
use unrestricted Hartree-Fock (HF) combined with an in- 
homogeneous generalized Random Phase Approximation 
(RPA) for linear lattice fluctuations jl^l in a supercell of 
size Nx X Ny with periodic boundary conditions. In this 
model, doped holes tend to segregate into the stripes due 
to competition between the magneto-elastic interaction 
that favors global electronic phase separation and elec- 
tronic kinetic energy that favors uniform carrier density. 
We note that other competing interactions can produce 
stripes, clumps and other inhomogeneities [3ll l3^ : the lo- 
cal, coupled charge-spin-lattice dynamics governing edge 
modes on these templates then can still be modelled with 
a similar Hamiltonian and RPA analysis. 

The output of the calculation here is the inhomoge- 
neous HF groundstate and the phonon eigenfrequencies 
and eigenvectors. From the phonon eigenmodes, we cal- 
culate the corresponding neutron scattering cross section 



J II, 



(2) 



where R;(i) — R" -i- d; -f ui{t) is the position of the 
l-th oxygen atom expressed in terms of the location of 
the unit cell origin Rj^, position within the unit cell d;, 
and time-dependent vibrational component u;(i). For 
phonon modes with u; [t) oriented along the correspond- 
ing metal-oxygen bonds, on the O^: sublattice d; = |a: 
and u; = xix, and on the Oy sublattice d; = and 
U/ = Diy. The scalar displacements can now be expressed 
in terms of the normal modes Zn as xi{t) = ^„ Oixi,nZn{i) 
and yi(t) = J^n'^yi.nZnit)- Making first order expansion 
in the oxygen displacements, we obtain 



(3) 



X- [(1 + nB)6{uJ - LOn) + UbS^UJ + LOn)]. 

Zraujri 



Here, = Z); e 'a^:,,™, and ub = (e' 
is the thermal population of the phonon mode n. This 
is a generalization of the usual neutron scattering inten- 
sity expression |33j for the case of phonons with a larger 
real space unit cell. We plot S'(k, ti;)/|kp for k-directions 
sampling longitudinal modes, consistent with the com- 
mon experimental convention. 

We first analyze the signatures of the edge mode by 
comparing computed phonon spectra for doped and un- 
doped cuprates (Fig. QJ. The predicted neutron scat- 
tering intensities are given in the original Brillouin zone, 



symmetrized with respect to the two possible orientations 
of the stripe. 
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FIG. 1: (a) Neutron scattering spectra and the ground state 
structures for the undoped (left panel in (a), and (b)), and 
hole-doped (20%) case (right panel in (a), and (c)). The neu- 
tron scattering spectra are symmetrized with respect to di- 
rections parallel and perpendicular to the stripe. In the lower 
plots, black circles represent Cu sites and red circles represent 
O. The radius of the circle is proportional to the correspond- 
ing site hole density. Arrows centered on circles show the 
magnitude and direction of spin. The green lines originat- 
ing from O sites indicate the magnitude and direction of the 
equilibrium O displacements in the presence of a stripe. 



The ground state configuration in the undoped case is a 
uniform antiferromagnet (AF), while, in the doped (20%) 
case, the holes partially segregate into a Cu-centered ver- 
tical stripe Is^l • The doping produces two low- frequency 
branches that split off from the band of extended phonon 
modes characterizing the undoped state. In Fig.[21we plot 
representative modes from these branches. Clearly, these 
are modes localized around the stripe. The dispersive 
lower frequency band centered around ^ 30 meV cor- 
responds to the oxygen displacements along the stripe 
and the narrow high frequency one at ~ 70 meV to 
the vibrations perpendicular to the stripe. (The ex- 
ample shown in Fig. right, is a localized asymmet- 
ric breathing mod e.l7...18.l .) Experimental evidence from 
neutron scattering [13,113 shows mode softening of about 
15 percent with doping in the cuprates, consistent with 
the softening obtained here for the perpendicular edge 
mode. Recent neutron measurements in twinned and de- 
twinned YBa2Cu306+5 ^3 have attempted to resolve 
the phonon spectrum both parallel and perpendicular to 
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the purported stripe direction[3^, but conclusive identi- 
fication has been hampered by the degeneracy and mix- 
ing of many modes near the zone boundary. There has 
been no search yet for the parallel edge modes which are 
predicted here to have a huge softening of ~ 60 percent. 
However, there are distinct anomahes in the longitudinal 
phonon spectrum of La2-a;Sra;Cu04 around 30 mcV that 
have not been studied in detail[l8j|. 

Despite the close agreement between the energy scales 
of the predicted and observed modes, it is important 
to note some deviations from the reported experimental 
data for the mode intensities. We speculate that these 
discrepancies may be due to the detailed local oxygen en- 
vironments: Either (a) more complex intra-unit cell dis- 
tortions, i.e. local oxygen polarizations (e.g., from local 
confinement or pairing effects) , which would require more 
detailed modelling than the constrained O displacements 
used here, or (b) more complex stripe configurations than 
the regular linear stripes we have considered. (Such tex- 
tured, e.g., checkerboard, patterns are indeed found in 
certain models 13^'). 
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FIG. 2: Localized vibrational eigenmodes for the case of a 
20 % doped Cu02 plane. The stripe is centered along the 
middle vertical Cu row. There are two localized branches: 
One corresponds to the oxygen vibration parallel to the stripe 
(low frequency, E = 14.8 meV, left), and one corresponds to 
the oxygen displacements perpendicular to the stripe (high 
frequency, E = 68.5 meV, right). 



We now investigate the structure of corresponding elec- 
tronic excitations in the charge and spin channels. We 
again use the inhomogeneous generalized RPA; however, 
due to the increased computational complexity, we must 
use smaller system sizes: A 3 x 3 system for diago- 
nal stripes and a 3 x 4 system for vertical stripes. In 
both cases the hole concentration is i. We have ana- 
lyzed the spectra, as well as the real-space structure of 
these excitations. For example, the excitation spectrum 
in the longitudinal spin channel is given by I{uj,q) = 
J2nMSziq)\n)\^S{Lo - En + Eo), where |0) and \n) are 
the Hartree-Fock ground state and an RPA excited state, 
respectively. The real-space structure of the excitation 
\n) is described by, e.g., {0\Sz{i)\n). This quantity shows 
how the expectation value of the z-component of spin 
on site i changes due to to the excitation of the RPA 
state \n) when the system is driven by an external field 
of frequency £'„ — Eq. 



The phonon spectrum for the vertical 3x4 case is es- 
sentially the same as in the 4x5 case in Fig.Q] with two 
localized branches near 30 meV and 70 meV. This insen- 
sitivity to doping is a characteristic of the edge modes. 
For the 3x3 diagonal stripe case, the phonon spectrum 
has only the high-energy edge mode at 70 meV, but no 
excitation near 30 meV. Note that the low-energy vibra- 
tional mode can only exist for a vertical stripe because 
the oxygen-mediated charge transfer along the stripe is 
suppressed in the diagonal case. 
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FIG. 3: 3x4 system with a vertical stripe. (a) Spec- 
tral weights of longitudinal magnetic excitation (solid lines), 
transverse magnetic excitation (dotted line) and charge ex- 
citation (dashed lines); blue lines are for (0,0) and red lines 
for (|7r,7r) momentum transfers. The scale for the transverse 
spin mode is 4,000 times larger than the one for others. The 
magnified spectrum for the charge excitation at (|7r,7r) is 
shown in the right panel, (b) The ground state, (c) Spin and 
charge configurations for the excited state ~ 40 meV; this 
frequency corresponds to the oxygen vibration parallel to the 
stripe. Here Sx = 0. (d) Sx configuration for the excited state 
~ 50 meV, where Sz = 



The spectra for magnetic and charge excitations for 
the vertical stripe case are presented in Fig.|3| The spec- 
tral weight of the charge excitation aX q — (|7r,7r) has 
a peak near ~ 40 meV, which corresponds to the fre- 
quency of the "half-breathing" oxygen vibration along 
the stripe (same low-energy branch as Fig. but with 
qy — it). This charge excitation produces dynamic dimer- 
ization along the stripe, which couples to the anti-phase 
oxygen vibration parallel to the stripe. On the other 
hand, the excitation near 70 meV is associated with the 
anti-phase oxygen vibration perpendicular to the stripe 
(same split-off branch as Fig.l^j), and corresponds to the 
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charge excitation at the Cu sites next to the stripe. There 
is also a magnetic longitudinal excitation at momentum 
(0,0) near 70 meV, which is coupled to the anti-phase 
oxygen vibration perpendicular to the stripe. The trans- 
verse spin channel mode at (|7r,7r) shows not only the 
Goldstone zero-mode but also a higher energy excitation 
^ 50 meV. From Fig.|31(d), this higher energy excitation 
likely corresponds to the spin wave branch folded back 
at the zone boundary (larger systems sizes are necessary 
to be definitive.) From continuity of spin wave excita- 
tions, it follows that there is an intensity somewhere be- 
low 50 meV at (tt, tt). It is possible that this mode has 
a relationship with the "41 meV" magnetic resonance 
mode[33, 113 ■ Note that there is no coupHng of mag- 
netic excitations with the low-energy vibrational mode. 
This is because in the Cu-site centered stripe oxygen vi- 
bration only transfers charge between the sites with no 
spin polarization, and hence does not couple to the spin 
channel. 
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FIG. 4: 3 X 3 system with diagonal stripe, (a) Spectral weights 
of longitudinal magnetic excitations (solid lines) and charge 
excitations (dashed lines); blue lines are for the spectra at 
momentum (0,0) and red lines for the ones at (|7r, |7r). (b) 
The phonon mode at 64.3 meV. (c) The ground state and (d) 
the excited state ~ 65 meV. 



Fig. ^ (a) shows the spectral weights for excitations 
in the diagonal stripe case. Here, there are several RPA 
excited states near 65 meV, which lead to a peak in the 
longitudinal spin excitation spectrum. The excited state 
that has the lowest energy of the states ^ 65 meV is 
plotted in Fig. 0](d); the vibrational mode correspond- 
ing to this excitation is shown in Fig. 21(b). As discussed 
earlier, unlike for the vertical stripe case there is no local- 
ized vibrational mode around 40 meV, and accordingly 
there are no RPA spin or charge excitations in this energy 
range. The only excitations of the longitudinal magnetic 
and charge channels appear in the 65 ~ 90 meV range. 
The transverse magnetic excitation shows only the Gold- 



stone mode. 

We conclude that the formation of split-off local 
modes in spin, charge and lattice channels, familiar 
in the polaronic physics of conju gat ed polymers and 
quasi-lD charge transfer salts|23, |2y, |33, is also ex- 
pected to be prominently manifested in the case of other 
doped strongly correlated materials, includingnickelates, 
cuprates, and bismuthatesQ, S H S H E 13 ■ The 
intensities of these signatures should increase with dop- 
ing (and therefore density of interfaces). In addition to 
ARPES and inelastic magnetic and lattice neutron scat- 
tering, experimental techniques that should be applied 
to study the edge modes include dielectric constants, and 
infra-red and Raman spectroscopy^^; resonant Raman 
scattering would be particularly valuable even for low 
doping levels H^. 
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